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ABSTRACT 


The occurrence of independent events at random in the 
plane, i.e. the formation of a planar point process, iS 
discussed. Both homogeneous anc nonhomogeneous processes 
are considered. A specific functional form for the parameter 
in a nonhomogeneous planar Poisson process is used to 
illustrate the development of test and parameter estimation 
techniques. The M ien finds application in the description 
Per biological phenomena as well as in search and detection 


problems. 
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i. INTRODUCTION 


Many probiems arising naturally in a physical sense are 
Seem sO comple chat tne Toentification and description of 
underlying mechanisms must use the tools of probability and 
ENGUNSUICS. oOMemor tne reasons leading to the requirement 
of using these tools are: 

Gai) the data base may be so large or complex as to 
Preeimde Identification of any driving mechanism 
without recourse to statistical analysis; 

(ii) if identifiable, the mechanisms may be inherently 
pussspDISUclc; or 

(iii) if identifiable and deterministic, the governing 
law which the mechanisms obey may be unknown. 

This paper is concerned with the use of statistics in 

mic UWGctitealweaLlon and mathematical descrintion of the spe ial 
distribution of events (occurrences). Included is the det c- 
tion and estimation of parameters which influence the 
desersiorstonmor this distribution. 

The area of concern here is a departure from those sta- 
tistical methods which have been developed to detect the 
effect of varying a controlled segment of the underlying 
mechanism. Among those methods would be the design of exper- 
iments, regression analysis, time series analysis, and 
analysis of variance. One goal of such analysis is to 
hopefully vredict the advisibility of pursuing some course 


Brsacetion. 





In the basic model of this paper, events are considered 
CORO CCUM ROSSON distribution in the plane. This 
Hearne navuralemede! ror the expression that ‘points are 
distributed at random'," [Fisher, 1972, p. 141]. The bi- 
variate Poisson process will be defined and then developed 
through the use of partial C erena e ace equations, 
a widely repeated procedure in the univariate case but 
neglected in the bivariate case. 

Initially a homogeneous Poisson process will be assumed 
to control the underlying mechanisms. Then trends will be 
introduced by defining the Poisson parameter in such a way 
as to make it be spatially dependent. This will be the basis 
for the definition of the non-homogeneous Poisson Process: 
Time inhomogeneity will not be considered. Thus, the data 
are assumed to be taken concurrently, i.e., the period of 
observation is short compared to any period of change of 
the parameters. 

Tests will be developed to distinguish betwsen homogeneity 
and non-homogeneity and the method of maximum likelihood 
will be used to develop estimates of the parameter in the 
homogeneous case and parameters in the non-homogeneous case. 
In the latter case, conditional likelihood techniques will 
be utilized to develop tests and estimates. Throughout, 
testing and estimation procedures will be based on a single 
realization of the process which consists of the number of 


events observed and their spatial locations. 





The problem Sr concern finds application in the estima- 
tiono Che density of trees in a forest; here one might be 
concerned with estimating the potential yield of nes from 
a given forest area where inhomogeneities arisemaue TO soll, 
Kramer patternes, Topopraphyzand other physical reasons. 

Another application might be in naval search and detec- 
wen prooblems. Bor example, one might be searching for a 
merchant ship in distress whose location is not known ener 
due to failure of the ship's communication equipment. Here 
the independence assumptions of the planar Poisson process 
may be valid, but not the assumption of homogeneity.  In- 
homogeneitics of location occur because of preferred sea 
lanes and physical characteristics of the ocean and 


atmosphere. 





ER BOMOGENEBOUS POTSSON PROCESS IN THE PLANE (HPPP) 


A. GENERAL DEVELOPMENT 

Comsuder a Ssuechasure process of events occurring in 
the plane (i.e., @ so-called point process) which is 
characterized by the assumptions 

T. There exists a finite positive constant A > 0. 

ge HO aheamnvemer= x - ib ana any sev Of non-overlapping 


1? , R, with areas Aj, ',A,, (in the 


usual geometric sense), the number of events occurring 


SS R 


in any region i, denoted N(R,), has a Poisson dis- 

tribution with parameter AA. which depends only on 

the area of the ao Aza and not its shape. Thus, 
n 

(AA) *exp(~AA, ) 


A A ¿O 


prob (N(R, ) n 
ieee, NCR, |, 1 = 15:29 "5k. are mutually indepe - 
dent in that NCR, ) is not affected by the occurren e 
of events in any other region or in any grouping cf 


the regions, G, as long as R flc = Y. Thus 


n, -AA 
k (AA, ) E i 
prob{N(R, )=n,, iz1,''",k]l- or ae (2) 


Definition i: If a process obeys the above assumptions it is 
called a homogeneous planar Poisson process (HPPP). 
For reasons of arbitrary shape the above basic definitions 


will suffice. However, under certain geometrical assumptions, an 





cqui va lenmdefinition Por the HEPP can be achieved in a man- 
ner similar to the development of the univariate Poisson 
PECceecmrnroughethesuse of partialed@tferentiajedifference 
equations. This is useful for the development of statisti- 
cal properties and will be very important in tke development 
of the non-homogeneous process. Such a development 2:50 
provides another phenomenological approach to the ho.zigenecus 
Poisson process, one which might arise through the struc- 
CU ea form instance. For illustrative purposes 
the following development will be accomplished using rectan- 
DIE Sans Note vhau the development is very dependent 
enscheseremerry involved; hence developments with OC nce 
geometries (e.g. circular regions) must proceed sc.ney7at 
differently. 

The underlying assumptions in the differentia: equation 

development will be 

ere exists a Linitespoesicive ccnstant& > 0. 

II'. For any region R* with incremental area AA, inde- 
pendent of one shape of the region except possibly 
as noted above | 
(a) prob {no event in R*) » 1 - AAA + o(AA), 

(b) prob (one event in R*) = AAA + o(AA), 


(c) prob (more than one event in R*} = o(AA), 


where "g(AA) is o(AA)" means iim ginn) EO sg 
/ AA>O 
specifically in rectangular regions the limit as Ax 
&(AA) 


BEL ver DOCHT EO TO Zero oí AXAyY is zero. 





III'. Tne occurrence 


of evints in R* is independent 


T 
Of the occurrence of events in any region R 


where RENR* = 
it will be shown cha? 


imoued by 1, IT amd III 


are equivalent and hence the i: 


>> 
> 
* 


Eise to a HPPP. Wearls 


ante. Also rfi impie 


(D prob INC R=? 


bs 


Terence definition of ol 


fi 
ps 
Lys 


(B) prob {N(R*) 


and (c) prob {N(R*) 


Iv 


The problem remaining in 


2} = 


2. 


aper 1 


and III' imply and are 


oe 


so that tre two sets of assumptions 


^ 
L 


cremental assumptions give 


Qu 


ft’ are the same,as are III 


ace by (1) 


1° 


2 
5 (AA) - 


l - AAA + 


2 


2 


SENE 
— AAXAYy + Ax Ay” - 


M AAA ar NOA), 


AA) given above. Also 


A AAO - AAA +...) 
= AAA eon Ay) 
" i -AAA 
E A = O(AA). 
ee 


order to demonstrate equivalence 


between the two sets of assumptions is to show that IiI' 


moles” Il. 


Consider a region R bounded by the co-ordinate axes and 


lines x = X* and y = Y*, 


with area X*Y*. Now extend the 





sides to x = X*+Ax and y = Y*+Ay (see Figure 1). Consider 
the probability of n events occurring in the extended 
region ' = PR ; 
eocuon. R R Ry + R, t R25 where: 
(a) R has area X*Y*, 
(b) Ry has area X*Ay, 
(c) R, has area Y*Ax, 
(a) R3 nas area AxAy; 
(a)-(d) imply R' has area X*Y* + X*Ay + Y*Ax + AxaAy. 


enemi ssumptions I!, II'Mand III' imply 


prob {no event in R, ? = 1 - AX*Ay + o(X*Ay), 
prob fone event in R, } = AX*Ay + o(X*Ay), (3) 
prob {more than one event in R,} = o(X* Ay); 


prob {no event in R^) - ] - AY*Ax + o(Y*Ax), 
prob fone event in Ro = AX) (4) 
prob {more than one event in Ry} - o(Y*Ax); 


and 


prob {no event in R3) = 1 - AAxAy + o(AxAy), 
prob {one event in R4) = AAxAy + o(AxAy), (5) 
prob {more than one event in Rh = o(AxAy). 


Moreover, statements (3), (4), and (5) are independent. 

It is noted that the above equations may have two 
desferent interpretations. For instance in (3), prob{one 
event in R,} - X*Ay * o(X*Ay) is interpreted to mean one even 


in a two-dimensional process with parameter A and area X*Ay. 


10 





jupe l. 





The incremental increase of a 
rectangular region. 


al 


A 


a. = 
i 





However, another interpretation would be to consider the one- 
dimensional(marginal) process of events projected onto the 
y-axis, in which case the parameter is AX* and tne incremental 
imverval has length Ay. 

For notational convenience, let P_(X*,Y*) denote ene 
probability that n events occur in a region with area X*Y*, 
iteeditiferentvial-ditference equations are written noting that 
n events may occur in an extended region by having n events 
in the ünextended region and no events in the extension, 

n-1 events in the unextended region and one event in the 


Sens ion, Ete. Hence 


Eta, rss ae) 


+ 


X vX) . ~ VX 
P _,(X*,Y%) + P)(Ax,Y*) 


re 


Pio (X*,¥#) + P(Ax,¥*) + 


= Pp (X*,YF)L1-AY*Ax] + Py (RAF) TAY FAX) + o(Y*Ax). 
(6) 


eumd anly. 
p (X*,Y*+Ay) = Pp, (00,14) [1-AX*Ay ] ES BEE an E 
+ o(X*Ay), ER 


and 


Te 





A yO = POUGCS,Y*)L1-AYTAx J[1-AX*Ay JE1- AAxAY ] 


Je 


Pq ARTE) LAY F Ax (1AXF Ay) (1-AAxAy) 


+ AX*Ay(1-AY*Ax)(1-AAxAy) 
+ AAxAy(1-AX*Ay)(1-1Y*Ax)] 


+ P FORSE Ay = AYA (1-dAx Ay ) 


N 


+ 1AX*Ay(1-AY*Ax)AAxAy 


+ (1-AX*Ay)AY*AxAAxAy] 


Cayó] (8) 


m 


¥ v#)[ yx 
Pa-3 X sY*) [A TX*Y*AX 


+ 


SENDER Er TR) 


Interpreting the above equations, the third term on the 
right hand side of (8), for example, states that there can 
be n events in the extended region R' if there are n-2 
events in R anc exactly cne event in each of any two of t 2 
added regions. That is, there can be two events in the 
added regions Ry Ry and R3 if one occurs in each of two 
regions and none occurs in the third region, i.e., one in 
Ry one in R, ande onec Res eve COllecting all termszert 
erder o(Axdy), o(X*#Ay) and o(Y*Ax), (8) reduces to 


= 





^ 


Pa Ax, Y*rAy)-P, (XF ,Y*) [I-AYSAx-AXFAy -AAXAy £A XEY*AXA;) 
+ P o (X*,Y*)[AY#Ax-2A X#Y#AxAy+AX*AY+AAXAY] 


r P QS Y*) DA XEY*AxAy H0 (AxAy ) £o (X*Ay )o (Y 8x) . 
Co") 


= PQOC,Y*)[1-AYYax]eP, y (X%, Y) [AY Fax ]+o(Y%Ax) 
+ POX* Y*)C1-axFay J+P | (X*,Y*) (AXP Ay ]t0(X* Ay ) 


- Py AE YE AP, (AE YE) AXAY FAP, (XF 1%) Axby 


2 


m OE RE SYS)WP. SX*,Y5) | AgRSSSO Gora 


Noting from equations (6) and (7) that the firet three 
terms on the right-hand side of the above equation ar: 
pP (X*+ax,Y*) and the next three are TO pet 
rewriting (8'), the result is 
E R a ME Wee). p (et M 
= AP (XE YS) AXAY*AP, y (X*%,Y%)AxAy (8") 
eyyy% * vk). ¥ V#)4D x vk ' 
T A XY [Pp (Xx ,Y*) eP _,(X BY JA A „Y* JAxAy 
+ o(AxAy). 


The definition of the second partial derivative with 


respect to two variables is 


14 


N 
1 


y: 





A Aim ÉAXtAxiythAy)-P(x,y*Ay) 


e Ay>0%%Y Ax>0 Ax 
= 145 Prax eR Gy), 
Ax >0 Ox 


uoe, tCranspesing the fist theee terms of equation (8") 


pm me riecnt hand side, Gee Gs, E 


PAV and enen taking 


mee double limit resu?2t- " 


2 x ve 
9 PROUGYT) LP (RR tar E, YE) AER [Pp (O, Y) 
roce n n=; 2 n 3 
IXOYy j u 
=2P EEE Ne] (9) 
imemsouutton to (9) (a oina cicferential-difference 
equation) is easily shew =. =. 


Baer) = ECO) EA. tl LO) 


m 
Cr. 
Te 


where K is an arbitrary constant. Special considerations 
pecmrcedcamnuor n - O,1 Since for these cases some of the 
terms -in (8") and (9) are not defined. Rewriting (8") and 


(9) while concurrently eliminating the proper terms leads to 


K(Ax*v*)l'exp(-AXSYX) 
n! 





Ne. EE So uo OO 


Since P (X*,Y*) is a probability statement and for any 
n 


vem reprom che number of events in that region must be 


some non-negative integer, the constant K is seen to be unity. 


15 





Hence (105) is equivalent to (1) which was to be proven. 

Thus tne two sets of assumptions imply the same things, 
nameiy that the number of events in a region has a Poisson 
Crs ue oUU LON E parameter proportional to the area of the 
region and independent of its shape and the number and 
position of events outside the region. Note that the 
formulation ezcludes multiple events, i.e., the occurrence 

DE WO Or more events at any polnt or on any line in a singie 


eaded region such as R in Figure 1. 


il 

Msc a similar derivation will go through for circular 
regions using polar co-ordinates, but there are differences 
in the special properties of the Poisson process as defined 


through assumptions I, II and III in differently shaped 


regions. These are discussed below. The differences in 


eicmowecilal properties of the mom-homogenecous planar Poisson 
Brocessessas whey vary with @iiterent geometmsdes are an 
essential element of the analysis of points (events) in 


the plane. 


B. TESTING DATA FOR HOMOGENEOUS TELANAR POTSSON PROCESS LHPPP) 
Given the occurrence and spatial location of n events in 

a rectangular region of area X*Y*, consider the problem of 

determining whether or not these points occur as realizations 

of the HPPP. Miles [1970,p.89] has stated a consequence of 

deninition l as 

Corollary. Assume a rectangular region R, with area Aj: 

Given NCR, ) = n and 0 <A, < o, the n points are independently 


i 
sodmunPtobpmiy distributed in Rie 


16 





Proof: Let A, - X*Y* where R, is a rectangular region 


al 
bounded by the coordinate axes, x = X* anc y = Y*. Label 
the n given points in any convenient manner. :.¢.. on the 
th 


magnitude of the y-component. Let (X,y) (1) denote the i 
labelled event. Consider an incremental region with area 
dxdy which has the property: prob fexactly one event in the 
incremental region of area dxdy} = P, (dx, dy) = Adxdy + o(dxdy). 
Consider now n incremental recangles dx,dy, ; ds ac ae 
placed in Ra. Ignoring probabilities of o(dxidyi), assump- 
Bons, IIT and III imply that tbe joint probability that 


the zen 


events 1n the incremenbal rectangle, ax, dy,» 
Mole *.ns and exactly n eventspoccur altogether in X*Y* 
is given by 
LAXXVE 
Adx, dy, ...Adx dy, exp( AX¥Y*}, 
Heswauime in terms of the density function, 


HAY py AN y 3) = Alexp{-AxtyY*}, 


mere £l...} 4s the joint density of (x,3) (4: los 1... a 
End the probability that the number of events in X*Y*ois n. 
The exponential term in the above expressions is an approx- 
imation to expl-(AX*Y* - i Aádx.dyi)l, i.e., represents 
the probability of no B within the region X*Y* but 
outside the incremental regions containing each event. 


By conditioning on the occurrence of n events in th 


region which are distributed Poisson with parameter AX*Y*, 


T 





\exp(-AX#Y*) 


Bee. (xy nA) =s u, n»1 
(1) (n) l (x*Y*)lexp(-AX*Y*) ` 
l n! 
n! 
Ricco —— (11) 
(X*Y*) 


which is the joint distribution for n bivariate uniform 
random variables ordered on one of the random variables as 
is shown in Appendix A. Note also the independence of the 
conditioned density from the parameter A, i.e. the random 
pusopwMe Ns a sufficient statistic for A. 

As a consequence of the above corollary, it is apparent 
that if the points of the HPPP, conditioned on the number 
Sie Nuts Observed to occur, are in fact ordered with 
respect to the increasing magnitude of the y-component, then 
no "information" is available about the ERE of the x- 
components, 1.€e., each of the n! orderings that can be 
induced on the x's by the orderings on the y's has probability 
1/n!. This is readily apparent since in the bivariate uniform 
case the two components were independently selected. Hence, 
Af (x,y) (y is determined by (AY y? fe. the points are 


labelled by the ordered y-component, then 


EIL 
Cs 
I 
> 
w 
NO 
v 
w 
D 


prob{X, = Xy? ES 


is Ehe on Xain- magnitudes ang 


where Xe k 


j) 


prob{X, = TE JOD SUD X57X (ty) » ke 1 dub] oom 


G XQ 7 Xa) n zT (12) 
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Hence if the x-components of the points ordered on the 
y-components exhibit any natural ordering then the x- and 
y-components nave not been independently selected and the 
apeielnvea process cannot be a HPPP. Fils will be the basis 
Om menye of the tests for a HPPP agatiist. S non-homogeneous 


planar Poisson process to be discussed later. 


Lemma: If the biveriate process is Poisson and the regions 
ES rectangular, then the projections of the events onto 


the coordinate axes may be shown to be univariate Poisson. 


Proof: Comsider the oecurrence of events imea rectangular 
region of area X*Y*. Then by III the occurrence of an event 
in an incrementel strip Ys independent of all occurrences 
emizside the strip. Hence the projections onto the coordinate 


axes give ríse to independent counts along the axes. 


_ (AYEx)Pexp(-1Y*x) 


PS O. 
E 3 7 n! 3 O < x < x5 
and 
4 gl - 
| = (A) exp. A) = 
P. (X#,y) = (AX y) ex (-AX*y) — n2 0,1, 
aa AS E 


which gives the univariate Poisson distributions with 
parameters AY* and 1X* respectively. 

Note here the inherent dependence on the shape of the 
assumed regions. In using rectangular regions equal lengths 
in the marginals reflect equal areas in the bivariate 


ais treDut ion. 
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If the regions were circular then vertical projections 
onto tne axes would represent decreasing area as the dis- 
tance from the origin increased. Since the occurrence of 
events is assumed to be proportional to the area projected, 
an actual HPPP would induce a non-homogeneous process on 
the marginals due to the distortion in the mapping. For 
clarification, refer to Figure 2. However, if the re ens 


Sue Circular chen radial projections could be MW. 1132048 050 


that the event occurring at (x ) in Figure 2 is repre 


99 
sented in the x-marginal by an event at X5. To define 
equal area projections in this case the transformation 
xo x* = x' is made, in which case a unit increase in x' 
derines the addition of a unit amount of area to the wa sicc-. 
For example, 1f a unit area is generated by a circle 
radius r = 1, then the area enclosed in the ring cr 
i <r SUS ne UNIDAS Sao ple femea in the wine 

of 2 as 5, Eco lmn general, yn <r a ie Get nes ve) 
polar coordinates a ring with unt area. 

Returning to the assumption"ot rectangular reglors, 
Ince characteristics of the HPPP are now available which 
can be used as the basis for testing a sample for belonging 
to the HPPP description of events in a rectangular region EH: 

(A) Independence of the x-ordering from the ordering 
on the y-componencts. 


(B) Univariate HPP (homogeneous Poisson process) in 


the x-marginal and, conditionally on n events in R; a 


20 





— 





Laure 2. Verri koa iena radial projections ofan 
event to form the marginal process. Shaded regions 
represent the deviations of projected areas arising 
rem. the reetangular proñeetion of seireular. areas. 
Thus, the shaded regions indicate the degree of 


non-homogeneity induced by the mapping. 


an 





enorm distribution of the distances to events. 

(C) Univariate HPP in the y-marginal and, conditionally 
Soma CVCives tt Ro, dWmeiomm distribution of the distances 
to events. 

Property (A) can be tested against general alternatives 
using a rank correlation procedure (or Spearman's correla- 
tion, see Pearson and Hartley [1966, Table 44]). Properties 
(B) and (C) can be tested by standard univariate methods 

as in Cox and Lewis [1966]. 

Note that in the above discussion the interest lies in 
the nature of the process rather than in specifically 
Gesewrping tne process. sede re rmimnation*tet te 
wemamecer A of the Poissen preeesses nota current objec- 
tive and it can be considered to be a "nuisance" parameter. 
Hence the conditioning argument above and the resulting 
Independence of the tests from the value of the parameter 
are justifiable. 

Nev te emo De Tee proband many ot 2 Type T error gener 


A 


Ancares tana for randomness, as be the corresponding 


B 
probably in testing for HPPP in the x-marginal, and a 
likewise for the y-marginal. Then the probability of not 
Suselkzreseerimertche HrPPenypothesis due to the randomness 


best ys l - o etc. Hence the combined probability of not 


A? 
falsely rejecting EPPP is 1 - prob {type I error} or 


arch) = (= ay) (1 2 a4) (1 = ac). 
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Therefore 


Ber 1 — (1 -= &,) (1 - ag) (1 = a) (13) 


is the probability of falsely rejecting a HPPP hypothesis. 
If through physical considerations one of the tests seems 
more or less significant than the others, the analyst can 
choose the weightings to so reflect the physical properties. 
Otherwise the values (and thus the tests) can be weighted 
equally. This need for the determination of weightings is 
ehewinherent disadvantage of a multi-level test. 
The individual tests proposed above will be briefly 
described. For the rank correlation test, consider each 
Xi from y) (ay which is ordered on the y-component. Also 
consider the ordered realizations along the x axis, where 


X, = X(4)* Then 


i=l 


S 7 n(n 


me 
6 E (i- (3),) 


e > (14) 


I 


where “DR is the position of x, in the x-ordered sequence, 


al 
is the rank Eco SUN S CNN 
The er Gistreours One ror P, can be approximated by 
fitting a distribution to its moments as discussed by 
Kendall and Stuart [1951, p.477]. The exact distribution 
of Po iS Cabulated in Biometrika Tables For Statisticians 
[1966, Table 44, p.23] for observed values of n between 4 


Zus aneeerne Introduction to these Cables gives approx- 


Bons for 10 <n < 20 and for n > 20. For 10 < n < 20, 
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P, can be treated as a product-moment correlation coefficient 
between normaviy distributed random variables. For n » 20, 
Pe Ie asiemmed” to be unit normal. 

Baar Sei megcne mazeimal distribution for HPP, two 
opa ace Cestorare proposca.srirst, the uniform conditional 
S Ued OR est aLa ns trends in the data. Ws stated 
in Cox and Lewis [1966, p. 153], "If the series has been 
observed for a fixed time t. (length X*) and n events occur 
m (0,t )((o,X*)), then the uniform conditional test is 
based on the variables U) - T,/t, {= X es) /X*}li=l,...n) 


conditionally on N being equal to n." The {brackets} are 


E 
Ec 
supplied to relate the material in Cox and Lewis [1966] to 


fas specific problem, and N- = n means the number of 
O 


posunsenceswob5erved is n, Nove that in the conditioning ~ 
of the realizations the "nuisance" parameters AX* and AY* 
easereliminated. | 

boron azvest based on the ordered inter-evenv spac- 
es is Used Co test Poisson apsenst steblonabry event 
processes which may be non-Poisson. For this test, Durbin's 
modifications of the uniform conditional test is used Cox 
& Lewis, 1966, p. 155]. Referring to Figure 3, Durbin's 
medication describes a transformation from the random 
variable X to the random variable T and then to the random 
variable S. 

Let T 


ze 


describe the "times to events" in a Poisson process, then 


Eger 1 = Me... Ms 


the Tr MEL mlaresudependent exponentially 
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Breunre 3, The generation of the transformed variables 


Si from the original process Xs): 
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distributed random variables with parameter A. If the 
T,'s Encaccnenwordered cand the S158 are generated as shown, 
then the S'S are independent exponential random variables, 
where S; has the expectation 1/((n+2-i)1). Also the trans- 
formation A = Gen defines independent identically 
distributed exponential random variables with parameter 
A, and therefore ds = : J a ms det ines the 
times to events in a u process with parameter A, 
and U. = ES, is the statistic upon which a new 
uniform en Bestes ed 

POS sc sesion Laden tala ed das Che Uniform conditional 
test is more powerful when testing for trends while Durbin's 
modification is relatively more powerful in testing against 
stationary event process alternatives. However, these 
tests are not independent of each other and thus cannot 
be combined as in (13). 

ancla os ne labo ve procedure, the region 
of concern may be partitioned into several sub-regions and 
the number of evénts in each subregion used as a basis for 
E tesina. This method is discussed by Kendall and Stuart 
[1951, pp. 574-5] who mention the problem of choosing the 
acne" partition, adding “Whether a particular partition 
has statistical interest depends on the purpose of the 
analysis". Due to the underlying uniformity of the condi- 
lana distribution, this problem reduces to the selection 


of the number of regions which are then used to form equal 


Tea SubD-regions. 
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Another alternative to the above testing procedure is 
the evaluation of the sample product-moment correlation 
C eent under the bivariate uniform distribution. The 
procedure is discussed by Kowalski [1972], but unfortunately 
his discussion does not address the bivariate uniform 
distribution. Kowalski makes two points very strongly: 
wutestly, the distribution o$ (the sample MU 
correlation coefficient unce: nen-normal assumptions) may 
differ from its normal-theory form and, secondly, we may be 
in a situation in which p is a poor measure of association." 
Hence, if the exact distribution for r under the bivariate 
uniform distribution were known, then an exact test for the 
HPPP (given the occurrence of n events) could be devised. 

Durbin [1970] has alse vroposed distance methods for 
testing bivariate distributions. he process herein described 
is well-suited to the methods Durbin uses since he first 
te nS orms the observations s0 that they occur uniformiy 
Smet ne Wait square. Hence the@navural trensformatior 
x' = xX/ "and y' = y/Y* avoids the problem of possible lack 
of uniqueness which is the central objection to the use of 
distance methods. These methods allow the analyst to adopt 
Durbin's bivariate analog of the Kolmogorov-Smirnov tests. 
The advantage of this method is the elimination of diffi- 
culties concerning multi-level tests and partitioning 


tests. 
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UDeEbests desewEbed in thisfsectilon are very general in 


Nauvure, ILe., 


ice prOceSS teh er P 1s vested acainse 


Hy: Cnenprocess is not HPPP. 


Rence the alternatives being tested against are multitudinous. 
Mamet 25 desired to test a realization as being from a HPPP 
Fan. a specific form of departure from HPPP, better tests 
may be defined based on the nature of the specific alterna- 
tive. For instance, one such departure could be non-homoge- 
Bun  ı.e., where A is not considered to be constant but 
ae Me fünction-of location this subject- -is considered 

in chapter III. Another departure might be in the nature 

oí the process itself. For example, events may occur 
according to some fixed plan in which case the process is 
deterministic and thus non-Poisson. A process may develop 
arena the occurrence of an event prohibits the occurre ce 
oí another event for some interval about itself, in which 
case events are not independent of other events and are 

thus non-Poisson. 

LG must be remembered, however, that tests against 
specific alternatives may ignore some features that a more 
general test would detect and thus each individual specific 
test applies only to the specific form of departure being 
considered. 

Moreover, in all reasonable stationary alternatives, 


it does not seem possible to derive the likelihood function 
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of the observations. One thus cannot derive exact tests. 
For tests against specific alternatives based on distance 
methods, see Holgate [1972]. Tests based on spectra are 


discussed by Bartlett [1964]. 


DS IMULATING A HPPP 

Suppose one were concerned with searching for submarines 
which are assumed to be dispersed in such a manner that the 
locations at any moment are generated by a HPPP. If one 
Search procedure is to be selected from many proposed search 
procedures, then a possible manner of comparing tne effec- 
tiveness of the proposed procedures is to utilize each pro- 
Sere against several simulated dispersions. In such a 
simulation, the only "variable" which would be of interest 
Mete che procedures, so all variables seen as @etTecction 
amare lass1itfication parameters, facilities available, etc., 
would remain constant. Another problem which might be 
considered would be the effect of the change of such param- 
er Bons search procedure selected (l.e., a sensitivity 
analysis of the procedure to assumed operating 
characteristics and facilities). 

Brzthesinitial remarks of Section B above and the 
statement of equation (12), several methods of artificially 
generating realizations of a HPPP can be determined. These 
methods may then be utilized to simulate the HPPP. 

Assume that the parameter AX*Y* is given. To select 


the number N of events to be observed in the region with 
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area X*Y*, generate a random number U distributed uniformly 


Ene cet N= nh if 


Mel. 3 De 
y (AX#Y#)Iexp(-ıXty®) HT. de (AX#Y#)Iexp(-AX#Y#) 
r 1! Sese 1! 


(15) 


The summations can be evaluated using either E or Gamma 
Integral Tables [Cox and Lewis, 1966, p.24]: 
n-1 i -u 
> ne 


uM 


2 
A RISOD Dos Saas) 


o n-l -v 


V e 
Ta dv. 


CASI 
u 
Post Consider a random sartable X distributed uniformly 
over (0,X*), denoted X n U(0,X*), and another independent 
random variable Y n U(0,Y*). As realizations of each 
ee ncemsweriable are generabed, number them chronologically, 
mae. order of appearance. Generating n (as determined 
above) such realizations of each random variable yields 2n 
numbers: Xj eeta XoY] «Y: 
The final problem remaining is to select a scheme for 
mating the x- and y- realizations to form ordered pairs 
Dee CC onstitute the realization of the HPPP. A few 


such schemes are enumerated: 


n 
l. The sequence S Ormsa RPPP: 
n 
2. If the y, are ordered cov *9oPm won ae then the 


sequence Sx, Y 0932? forms a HPPP. 
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3. Similarly, <(X gy) 4)> forms a HPPP. 
Hw Additionally, any random permutation of the x, in 2, 
the Y4 in 3 or either random variable in 1 can be used 


Comi orm a HPPP. Thus <(x > forms a HPPP, etc. 


p 
The goal of the simulation and the purpose of simulating 
Wiese process aS a part of the overall analysis must now be 
Eouswdered. If during the simulation it is desired to 
generate independent realizations of the process, then each 
iteration must involve a selection of n, the drawing of 2n 
uniform variates and the mating of the variates through some 
scheme such as those outlined in steps 1-4 above. On the 
ser hand, if it is desired to utilize variance reduction 
techniques, then for any drawing of 2n random variates 
pEsehemes could be used for the mating process. Here 
independence is lost immediately and this loss must be 
balanced by some gain elsewhere in the analysis. 
D. ESTIMATION AND TESTING FOR THE PARAMETER FROM A HOMOGE zOUS 
BrRANAROPO PS SON "PROCESS (HPP2) 
Tf the hypothesis that the process is HPPP with some 
unknown value of the parameter A is accepted, one might 
like to Bras @ point estimate wer commaudesee interval 
Soma te. tor A. or to test that the process has some given 


parameter A Note that the parameter A, which was considered 


0° 
to be a nuisance parameter in the previous section where 
the structural aspects of the process per se were tested, 


Mewsspeesites the process completely. 
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Since, as Waszseeneinsseerjion Pat is possible to 
set up Te Jat Probability density function-of- the 
Sbservations in a HPPP, point estimation of A can be 
based On the method of maximum likelihood. Note, however, 
mac edach observation consists of a single "look" at (or 
realization of) the process rather than n observations of 
single random variable. Since it is a stochastic process 
the observations are not independent and identically dis- 
merpuved, Hence the usual justifications for maximum 
likelihood procedures are not valid; see Brown [1972] for 
extensions of maximum likelihood theory of estimation to 
ganizations oÍ a Poisson process. 

Using the results of Brown [1972], suppose that n HPPP 
events are observed to occur in a rectangular region of 


Sas Y. From (11), for n > 0, 


~1Y¥¥V¥F 
no AX*Y als) 


ti 


L = PU) y e Y) yy 3) A 


in la AX*Y*, (0 ST) 


If n £ 0, this function is -e at A = 0 and X = » and since 
p - = - X*Y*, the slope of the function decreases 
monotonically from e to -X*Y*. Thus In L has a unique 
maximum at the point where 3 = 0. Setting this 
derivative equal to zero yields a unique maximum likelihood 


point estimate for A as 


d= TEE» (n > 1) (25 
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— T E 
where À is unbiased (since E(A) = zur = \) and has variance 


A/X*Y*, Note that as the observed area X*Y* becomes large, 
the variance of the estimate becomes small; thus, by 


Chebyshev's Inequality [Lamperti, loeo p. 20] 





NEC 2S  —— ., (a»0) 
and as X*Y* ^ o, 
PLA - Al > a} + 0 


for all positive a and hence \ [(uopcestolunoprobabsdit;y. 
emat ter Statement is equivalent to the assertion that \ 

is a consistent estimator for A. Also since the variance 

of A UND VASA \ has an estimated variance X/ XA YA = n/(X*Y*)*, 
and an estimated standard error of fn/X*Y*. 

Ene above method 1smmot applicable. In this 
case, it might be preferable to give a confidence interval 
estimate for A. Specifically, a one-sided test alternative 
is used to generate a test for the assumed value Anull using 
as an acceptance region only n = 0. Intuitively, Aa will 
be small enough so that er < 1 (i.e., the expected 
number of observed events is less than 1). The nypothesis 


to be tested is Ha: à = A USA A Defining 


0 null 1 muli- 
a level of significance a from (16) by 


es 
SEIN OR SS, Teure : 
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the hypothesis Bo is accepted at the level a. Conversely, 


O any elven value of a, N dd may be determined by 


N Ko rM E 


mull 
X =~ Iln (l - a) 
null XRrY% A 
where the b thus determined is the largest value of A 


that tne test will accept at the a level, given that n= 0. 
enano the case of n> 1, in order to test that 


tne parameter of the process has some given value A assume 


02 
that n events from a HPPP are observed in a region of area 


X*Y*. The hypothesis to be tested is HL: A = A, against 


0' 0 


the two-sided alternative E A A ho although one-sided 


caleta vives can also be considered. Since N is a random 
variable taking on all nonnegative integer values with som 


positive probability for any A there is always some 


0? 
possibility of an observed value of the random variable N 
(the observation being denoted n) falling outside any finice 
range of values. Thus a region (n^ ,n*) pmnusawDccsspmeedgpPqed 
such that iP N lies in the region the hypothesis Ho is 
accepted; otherwise the hypothesis is rejected. The level 


one test is the probability, given A ="’., that N 


0 
- + 
falls outside the region (n ,n ). 
Since the test has been defined to be two-sided, the 


level is split into upper and lower levels a and a 
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4 = 
so that a =a +a. The procedure must consider values of 


À < Ag as well as values of \ > hos To proceed, it is 


necessary to define 


+ + + 
ee i NS (18) 
_ “dX X#Y# 
© (ApX#Y#)Je 0 
= à . A UUO 
j=n* J} 
and 
a = Pie ne, re (19) 
= ¥V# 
A (AgX*15)J exp (-^o* Y*) 
= p : 
j=0 if 


Thus, for a given SUN an n* naewwe"determdined such tnat the 
statement (18) just holds. Also, for a given oa , a n may 
be determined such that (19) just holds. 

The null hypothesis is accepted at the a level if the 
observed value of N falls between the two prescribed limits 


Bes n ), where prob{N £ (r^ Hl = a. Note that as 


(n 
stated, the result is indeterminate since a, once given, 
leads to many values for a and a@ =a - a” which satisfy 
the given a. The manner of selecting a* and a^ must be 
stated. Arbitrarily it may bedesirable to set a* Jo 07 2 


Reemupwocacally, this choice of a symmetric acceptance region 


is reasonable since as n increases, the distribution of N 
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is approaching the (symmetric) normal distribution. The 
choice of equal a” and o” may not be reasonable, however, 
for small AgXT YT since the™ Potsgon distribution is 
positively skewed. 


The statement probíN Z (n^,n*)|A = àn} = a is the result 


ro 


of the test of the hypothesis H à = A, at a given, fixed 


0 0 
fever og. It is this result from which one must usually 
draw conclusions regarding specification of the process. 


If the information thus available, i.e. H, is rejected 


0 
or accepted at the pre-determined a level, is deemed 
meuni cient for the purposes of a decision maker (for 
example) then another possibility is that the post-analysis 
information may be extended by determining for each obser- 
vellón the exact a, a at which the hypothesis would have 
been rejected. The decision maker is then left with the 
problem of the determination of his own level of significance, 
possibly based on his CUT grasp of the problem and 
its significance in a larger frame of reference.  Qnce he 
has determined his preferred significance level, the hypoth- 
esis is rejected OAC C Preden Ne SPEC Mied leve Mey 
comparison with Qe. Thus the decision maker has gained 
some influence over the analysis but has had to pay with 
some time to reflect on the problem at hand.  Alternativeiy, 
he can use d informally as am goodness of fit" of the 
hypothesis. 

Berne (18) and (19), the significance test is defined 


conventionally [Cox and Lewis, 1966, p. 30] to be: the 
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hypothesis A = Ay would be accepted at the level of signi- 
ficance a in a two-sided equi-tailed test if the observed 
Enbeucof events, n, is such that n, when used abternatively 
in (18) and (19) (i.e., is assumed to be one or the other 

of the end-points of the acceptance region), produces Ot, 


acme solution to 


Be) = 2min(P,(n;3g),P (n;A19)] UAR (20) 
Note that each observed value of n generates a new d vor 


any assumed Ag, hence a, = a (n, E For example, 


(e? 
MSOS20)—= 09596, P(20;20) = .7628 and P(15:20) = .1332. 
It can be seen that the fixed level procedure is 


computationally simpler, since for a specified a and A the 


0? 
interval (n^ ,n*) need only be computed once while in the 
latter procedure d must be recomputed following each 
observation of N. 

The inverse of the above approach which utilized the 
two-sided equi-tailed test of significance for a given 
value Ao leads Aco the determination ot "cent i dencesminGemia! 
estimates of A. Given that n events are observed, it is 
required to determine some limits on the range of A such that 
the true ae value A* lies within the stated limits 


ana probability 1 -— a. That is, it is required to 


ee X "OD and a A (N) such that 


PLATON) < A® < X(N) IN = n}=1- a. (930) 





Using P{N < n]\ = 4°} « 1 - a? to define a A' as the 
greatest A such that equality just holds and similarly 
using P{N < njA = 4} = a to define a A^ establishes the 
limits such that (21) holds. For a proof of this, see 
Brownlee [1965, p. 121]. Note that for each realization 
of N, a new ordered pair (A ,A') is defined so that the 
ordered pair is a function of a random variable and hence 
is itself a random interval. The procedure only states 
that for (1 - a) x 100% of the observations the true 
parameter 1% will lie within the limits selected. The 
limits for observed n from 0 to 50 are tabulated [Pearson 
and Hartley, 1966, Table 40]. 

Hor a normal approximation to the confidence interval, 
Cox and Lewis [1966, p. 31] define the upper a point of the 


mus enormalt dascPrbution as COD and give the relationship 


e x 
PE AXFY 


rob{-c 
: OGxsvs) 1? 


< C44! =l-a, 


la t 
2 


2 


the relationship being correct as AX*Y* > o, The confidence 
limits thus obtained are, to a second degree of approximation 


using a continuity correction and the estimate o(1) = Vn/XFYE, 


I 
met 5€ a t Cjan: 
2 


MP MW 


For example, if 50 events are observed from a HPPP, the 


exact .05 confidence interval is 37.11 « AX*Y* « 65.92 
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[Pearson and Hartley, 1966, Table 40] whereas the normal 
approximation gives 37.79 « AX*Y* « 66.07. 

¡Mere also exist N approximations Tomche syammnleance 
tests and confidence intervals [Cox and Lewis, 1966, p. 33; 


Brownlee, 1965, p. 173]. 





III. NON-HOMOGENEOUS PLANAR POISSON PROCESSES (NHPPP) 


Pee GENERAL DISCUSSION 

If the stochastic process described above is generalized 
desa Now the probabilistice structure of the event process 
tombe dependent on tne location of the events, a non- 
homogeneous planar process is evidenced. In the simplest 
such case a non-homogeneous planar Poisson process (NHPPP) 
Ases ia. an the definition of the Poisson process given 
aeoıe assumption I is modified to become 

I Rere exists a’Desitiye, finite function CRS 20305 
Also note that II is changed by the fact that the number of 
events in any region is not only a function of the area 
of the region, but also depends on the location of that 


region within the universe under consideration. Thus A 


is now expressed as À = A(x,y), and assumption II becomes 
TE., prob{N(R, ) zn) 
ACA, )}" exp{-A(A, )} 


n! 


where 0d |) ae f Alx,y) dxdy the symbol A^ implying the 


A 

integral over Ks area and A(x,y) is ea to be continuous 

over R, (with area Ai) so that the integral is valid. 
Assumption III remains unmodified, i.e. events occur 

independently of any other event or collection of events. 


Under the additional assumption that A(x,y) is continuous 


"Nin the region of “eonsideration, the incremental 


40 





development of Chapter II may be extended to achieve a 
description of the NHPPP. Additionally the continuity 
assumption on A and the definition of the parameter in the 
process as an integral over A eliminates re cere ley MES 

of line discontinuities, although there may be cases where 
this is an important component of the problem. This problem 
is not considered here. 

Bejenring back to Figure 1 in Section II-A, Consider 
specifically the incremental strip defining region Ry > IR 
the strip is divided into n sub-regions of equal area by 
taking n equal increments along the x direction each of 
lenesch 6x, then, under the assumptions on the behavior or 
A(x,y), the process in the ith sub-region can be approximated 
by a HPPP with parameter A'- Mxsy) 5 where (x,y). is an 
arbitrary point in the ith sub-region. Specifically (and 
arbitrarily) the Tower" left point 1s chosen for the 
succeeding discussion; thus thewparameter for the first 
sub-region has parameter à = A(0,Y*). Continuing, the 


probability statements for occurrence of events become 


n 


P, (x,y) A(0,Y*)6xAy * o(6x,Ay) MEE ES xe 


I 


P, (x,y) Nope wis o MELUO MA DEEEOSCE xe E292 


P, (x,y) Er or VA ON MEE Go ay) een 


where j = 0,1,...,n-1, Y* < y < Ay and ndx = xe. 


U 
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since the probability of more than one event in Ry 


is o(X*Ay) the probability statements above are additive 


and 


n-1l 
prob {one event in Ry) = E A(j6x,Y*)6xAy + o(X*Ay). 
4=0 


Bu une limit as n > ©, by the definition of an integral 


y% 
{ S A(x,Y*)àx)Ay + o(X*ay). 
0 


prob {one event in Ry 


By similar argument, 
yx 


{ f A(X*,y)dy) Ax * o(Y*Ax) 
0 


prob fone event in Ro) 


and 


prob {one event in R3} A(X* Y*)AxAy + O(AxAy). 


By comparison with equations (3), (4) ana (5) the above 
statements lead to definitions for average parameters for 
R5 and Re as 


each ofi the regions Ras 


ere") : x HE 

Dopo) E de FO way, Dez 
and 

Aa (X* ,Y*) = ee. 
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Using these average parameters, the equations (3), (4) 
and (5) are generalized, resulting in the following 


statements: 


prob [one event in R,) A iy + o(X*Ay), 


1 


prob {one event in R,) AS (X8 ,Y8)Y*Ax + o(Y*Ax), (23) 


and 


n 


prob {one event in R3) Az (AF, Y) Axby + o(AxAy). 


Using the result (22) as defining the parameter in each 
of the incremental areas in Figure 1, equations (6), (7) 


and (8") become 


Pea, Yt) 


n 


X vx a) VF X vX|YY vx 
P(x 3 ben: ÀY AxJ+P__,(X Sot JAY AX 


+ 


o(Y*Ax), (24) 


X yx š y% LA XX X vXYY YX, 
ECC PO eee ye OA TA y 


+ 


o(X*Ay), a) 


and 


P, (X*+Ax,Y*+Ay) pP (1%, ,Y*+Ay) + NS) AS 


a 


AZAxAy[P, (X*,Y%) - pen COR ] (26) 


+ XA A OE 2? (X#,Y*) 


m -1 


x + 
ige (D. CTS 


+ 


O(Y*Ax) + o(X*Ay) + o(AxAy). 
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Rearranging terms, dividing by AxAy and taking the 


double limit as Ax + 0 and Ay > 0 yields 


a°P (X*,Y*) 


FEN SRS ASCP (X*,Y*) Pa ea (27) 


4 Xy * X vyX). x vx Xx vx 
Ay AA Y BEZ N) 2P, Q AT. O ,Y*)] 
Dessen, together with the boundary condition that Ra Sra 


probability statement, gives 


P (X#,Y*) = Ce exp{-A(X*,Y*)}, n = be 
where 
Kay: 
AGEY*) = SS Alu,v) dudv. (29) 
; 0 0 

Thus, che number of cvonce occurring in a region 
Bounded by the coordinate axes; x = X* and y = Y? has a 
Poisson distribution with mean given by A(X*,Y*). The 
mean can be considered to reflect the cumulative effect of 
A(x,y) in the region of concern. 

If n events from à NHPPP are observed to occur in a 
rectangular region defined as usual with area X*Y*, and the 
events occur at Oy) ex)» l'z01. ns the labelling done 
Onschemmaenitwude of the y-combonent. then the joint density 
of the events and the probability that the number of events 


in X*Y* is n is given by 
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N ACOx y) (4 exp {-ACK* ,Y*) I. (30) 
j=l 
Note that (30) is a direct generalization of (16). 

Hence the NHPPP can be described in a fashion similar 
ger che HPPP, but the expressions have acquired increased 
complexity due to the necessity for the inclusion of 
integrals to define the parameters. The degree of added 
complexity is dependent upon the choice of the specific 
"funetional form for A(x,y). The next section develops the 


expressions for one specific form. 


DAA SPECIAL CASE 
To consider the location dependent type of process, 


a particular form for A(x,y) is chosen as 
Alx,y) = exp fat Bx + yy + öxy}. (31) 


Newe that if 8x * yy + óxy changes very little over the 


uneco m terest of x and y, then 

ASA = (1 + Bx + yy + öxy)expla}. (32) 
uer e lationships may be used, however, they may cause 
necessary and untidy restrictions on the values which the 


constants a, ß, y and $ may assume. In particular, A(x,y) 


must be greater than O and the bivariate exponential 
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polynomial (31) ensures this with no restrictions on the 
range of the parameters. 

Additionally, algebraic manipulation of the form reveals 
that the curves of A(x,y) » c, ec a constant, include some 


interesting properties. 


E o - 0, then In A(x, yo = € is a family of straight 
lines in the plane, intersecting the x-axis at an angle 

8 = tan+(-8/y). In this case a clock-wise rotation of 

the coordinate axes through an angle 8 would give an 
Expomentral function of y only. 

2. If 6 # 0, then 1n A(x,y) = e describes a system of 
contour lines which form a hyperbolic paraboloid with a 
saddlepoint at (-Y/6,-8/8) as is illustrated in Figure 4. 
itemay be helpful to interpret the Figure ir terms of a 
section of forest which has been sampled. The line r 
describes a possible direction of steepest ascent (DSA) 
whach passes through or near the region being sampled. 

Wis DSA may not be a topographic feature but rather a 
mathematical expression for a possible increase in density 
of trees along some line. Obviously, there may exist a 
SERON correlation between this mathematical DSA and some 
topographic features. Note here that along the DSA maximal 
values for A(x,y) are found in the sense that departing the 
DSA at right angles leads to decreased values for A(x,y), 


Ie. decreases in the forest density. 
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Figure 4. Contour lines for A(x,y) = expla+Bx+yy+8xy). 
NOS ROTOS Sana Cao De a descräpee 
(hatched.) Here B/y = 2 , 8/6 = Y and 


all coefficients are positive. 
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3. The exponential form can be extended with little con- 
ceptual difficulty, but possibly greatly increased mathe- 
matical difficulty, to describe a much wider range of 
possible circumstances. For instance, it is reasonable 
to assume that the DSA line will bend; hence terms such as 
ex“y and Exy^ and higher order may be included in the 
ponen. 

For the special form of (31), the cumulative or integrated 


intensity function A(X*,Y*) is given by equation (29) and 


becomes 
p. LAC). NEM n 8 
A(X*,Y*) = exp(o-8y/5) = O 


(33) 
- EL{(s+ove)t) + EMO, 


where Ei (-) is the exponential integral and Ei(-) = ce + in (:) 


T 


Tp » Where c = ‚577216 Ys a constant, as defined 1 


i=l 
Jahnke and Emde [1945, p. 2]. 
The likelihood function for the NHPPP may be develor :d 
In a manner similar to that used in the discussion of the 


HPPP. The discussion leading up to (16) is modified by the 


fact that the parameter is location dependent resulting in 


n 
L = exp{-A(X*,y*)} I 


WEN), (in > (34) 
: (4) Z 


T 


where (X59) (4) is a labelling of the coordinates of the n 


Pointe events. Thus, 
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ln L = -A(X*Y*) + 
al 


tus 


IL) ) 
; BR a 
follows directly. For the special case of A(x,y) given by 
(31), equation (34) becomes 

n n n 


ln L - -A(X*Y*)48 X x,ty © y.+6 X x 
ee a 


; ,y,*no (25) 


where Á(X*Y*) is given by (33). 

The above joint density, or likelihood, function pro- 
vides a functional form which may be manipulated to accom- 
Pishi the two principal concerns- -of the analysis of point 
pnoeesses: hypothesis testing and parameter estimation. The 
obvious null hypothesis is Hg: 8-2 y = 6 = 0, in which case 
the nonhomogeneous Poisson process 15 being tested for homo- 
geneity since a non-zero a yields a constant parameter 
à = exp(a) > 0. Should the above null hypothesis be rejected, 
— the analysis proceeds to develop estimates for the 
parameters 8, Y and 6. This phase of the analysis may 
procesa differently depending on now many and which of the 
parameters were tested as being different from zero. The 
complete, and most complicated, situation develops when all 
parameters are determined to be non-zero. Testing of 
parameters is the topic of Chapter IV while Chapter V 


discusses the estimation of parameters determíned to be 


non-zero as a result of the testing procedure. 


49 





IV. TESTING FOR NON-ZERO PARAMETERS 


A.  PRELIMINARIES 

it is desired to formulate a method for testing the 
data (i.e., the number of events and their locations) in 
order to determine which of the parameters in the model 
given by (35), specifically a, B, y and 6, are non-zero. 
Note that three assumptions are dnce at the outset: 
first, that the NHPPP model is valid; second, that the 
testing for homogeneity in Section II-B led to the rejec- 
tion of the hypothesis of homogeneity; and third, that the 
physical phenomena can be modelled by the NHPPP — by 
(31) with the parameter A(x,y) given by (31). 

Testing the Poisson hypothesis per Se when Thee eumiet Ton 
A(x,y) is not known is a compound problem which will not be 
considered here; It is analagous to the compound problem in 
regression analysis of testing both for an unknown regression 
Micron trend tor independent equate variance errors. 

From the third assumption, the likelihocd function for 


Wicmeavasts given by 
L = expl-A(X*,Y*))5% exp{Bix, + yZy, + SExy, } ESO ) 
where, for clarity in the future development, & = expfa}. 


Conadi tonilne on the occurrence oil n events, n > 1 and 


defining LE(X,y) (qyoe ees (XV) (ny l ns rCxsy) I 2n) leads to 
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! ex + e 
L NL exp(BZx, * yIy, * 6Xx.yi) 
prob (N=n ) Y X n 
( f f expíBu + yv + Suv) dudv) 
0 0 


L(n) = 


(37) 


where L(n) is read "the likelihood function conditioned on 
puemoecurprence of n events." Note that conditioning on 

the number of events observed has resulted in an expression 
which is independent of the parameter E(or a), i.e. for 
Pen Bs y, ana 6, nm is a sufficient statistic for a. This 
is convenient because a here is a "nuisance" parameter since 
Pies tems o interest are those which would indicate non- 
Mamo aeneloy rather than the establishment of the overall 
rate of occurrence. Thus by using the conditional likeli- 
hood a may be eliminated and the testing can proceed for 
non-zero 8, y and 6. In other words the value of a should 
not influence the testing for non-homogeneity parameters. 

mum -c0MEcertainiy no departure from homogeneity cou. 1 

be evicenced and hence this case is covered by HPPP; see 
IB boe mEHence the case of interest is n » l. 

Physically, the model (35) gives rise to a parameter 

surf$ace A(x,y) which has the properties: 

CINE on v Z 0; ó £ 0: in AgSsENs a hyperbolic 
paraboloid superimposed on 
amu ItedN)Iane. 1.9. c0me 
"warping of Lhe’ vtilved 
pioneers enidencede 


(b) 6 #0; Y XAO; 8 = 0: In A forms a plane, tilted 
with respect to the x=-y plane. 


(ey 2 0- 60 X 0: In A forms a hyperbolic 
paraboloid. 


ou 





Km = vue e ln A forms a plane parallel 
Demenz an”. ı,e, a 
HPPP is evidenced. 

There are a number of possibilities for testing: 


Coat testeof 


against 
E at least one of the parameters B, y, 6 X O 


is a test for non-homogeneity which is more specific than 
a: in Section II-B and is easily derived by likelihood 
tacio techniques. 

(b) The above test is not of great interest; generally the 
Specific non-zero parameter is desired rather than just that 
at least one of the three is non-zero. This leads to the 
suescudeme Nwselectimg these NSubseb;»5vasprobiem 
which is difficult and as yet is unresolved. 

(c) The simpler problem is to assume an ordering, i.e. that 
if B = y = 0, the process is homogeneous (6 is then assumed 
to be 0) and if B or y is non-zero but 6 = 0, then higher 
order terms are assumed to be zero. However, if the test 
indicates non-zero B or y this may be due to an aliasing 
effect because of a non-zero 6. If further testing of 

6 = 0 against 6 # 0 reveals 6 # 0, then it may well be that 
the true situation is 8 = y = 0 but 6 # 0. The procedure 


to be followed will not discriminate this case. 
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MMieesame aliasing effect occurs im testing of 6 = 0 
against 6 f£ 0 where B and E, are non-zero and it is desirable 
vo perform this test without the effects of the non-zero 
B and y. These are thus nuisance Paramnevers sas Vac che 
case with a in testing ß and y. For the present model (35), 
one can eliminate these parameters because it is seen from 
pi mer ponential form (360) that tor any o, (n, EX, > » D 
NN octuoft Sufficient statistics#for (a, 6, y). Thus 
6 = 0 is tested with some function of EX Ys given n, Lx. 
and Ly; This statistic has a distribution independent 
phe parameters a, B, Y. 

The reason for basing the conditional test on LXV; 
meo hachis is (Conditionally) a SUmricient statistic 


Hor 6. 


Bez SBESTRIE TESTS 

Assuming that some ordering exists on the parameters as 
discussed in possibility (c) above, tests are performed 
Helene surticient statistics (ma Exi, Ly}, EX; Y) to 
determine if any non-homogeneity is evidenced by the data 
(1.e., through the statistics). This testing is more 
specific in nature than the testing encountered in Section 
II-B above due to the selection of a particular model. 
The set of sufficient statistics arises from this choice of 


a specific model to use as an alternative to homogeneity. 


De 





tne testing will assume the following sequence: 
(1) Condition on n and set ô = 0. Test Hoi)? B=y=0 
against a B #0 or y #0. Note that it would not be 
informative to test either ß or y as a separate entity since 
in the formulation of the model ß and Y are unique only up 
to an angle of rotation. That is, testing of B and y 
Jointly amounts to the detection of any tilt in ln A(x,y) 
with respect to the x-y plane, regardless of the direction 
om Che tilt. Failure to reject Holi) leađs to the assumption 
of homogeneity due to the assumed ordering. 


(ii) Rejection of Ho) leads to testing of 


Ho(í19: $=0,-2%<8B<ow and -o<y<o 
against 
H4 (43)? ô ž 0; -2 < B <% and ->œ < y <œ, 


The test thus specifies y and B as nuisance parameters. 
SIS ENTS mne ceS Sar y t Ome ONTT LILOn on N, LX 
and ty. to eliminate the nuisance parameters. 

iio. cond. troning on mM ang sevueme o = UO leads co 


(By)? n! expiBZx. + YE yy) 


L(n) = — m . 
(exp{ßX} = (exp{yY} - wie 


From this it is seen that the statistics (Ex, ,2y,) are 


Keondieionallv) Jointly sufficient statistics for R and y. 
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Under ie: 


Ex, /n > N(X*/2,x*%/12n) and Ly,/n > N(Y*/2,Y**/12n) 


and the statistics are independent (see Section II-B). Hence 
the expression 
Ex,/n - X*#/2 2 Iy./n - Y#/2 2 
AA y o 0. 
Mey yon gy aun 
is asymptotically x5. Rejection or acceptance of Fo(i) 
is based on the adherence of the calculated value of this 


sum to the xs drstributionge nn eH. ucmccepped ir this 


0 
sum has sufficientiy small values. Acceptance of Ho (1)? 
Ux red earlier, leads to assumption of HPPP; refer to 
Chapter- II. 

Following the rejection of Hoi) it is necessary to 


proceed with testing of H As can be seen from an 


E 
examination of (37), the complexity of the exact distribu- 
tion following another conditioning argument (i.e. on 

n, LX; and Ey.) is puechdibubive oO Never cor large sample 
sizes the conditional distribution can be approximated from 
the fact that Ex, /n, Zy,/n and Ix,y,/n, conda tioned on n, 
A cata normally distributed for large n. Thus the 
asymptotic distribution of EX¿Y¿/n, given n, Ex¿/n and 


2y¿/n, can be found from normal theory multiple regression 


results. 


DID 





Under the assumption that B = y = 6 = O, the trivariate 
normal distribution which arises is characterized by a vector 
and a matrix. The vector (y) of expected values and the 


variance-covariance matrix (I) are given by 


| X/2 i D 
= 2 = 

i h 

XY /Y "xy | 
and 

X*/1?n 0 X*Y/2ln 

== 0 Y^/12n xy“ /2ln 
X*Y/2lin XY^/2Un . TX*Y^/1lln 


from which Pio = O and 014 * 054 = 0.65465. 


In the model given above, 


= Q; -0 < R < œ; -0% < y < 0 


aan) ° 
is to be tested against 


6 #0; -œ < B <œ; -œ < y<, 


aa)" 


SNE e EX: Y y ls a sufficient orario OA NENA 


and LY; are given, che test cankbe Dasced on LX.yi-. Its 


asymptotic (conditional) normal distribution has mean 


um ; . 
Hay and standard deviation Oxy given by 
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= o Hd E 7 
E Uy 5, JE (Ex, /n x/2) + 
“xy 
+ _ 
D P53 (Ly, /n Y/2) 
y 
and 
il 
i 2 DES 
xy ^ ?xy (o7 834 7 0234) 
Duy,/m = n 
E e AAA As distributed as a unit 
OLCI) Er 


normal variate and is accepted if this statistic has 


H 
0(ii) 


Surticientily small values. Failure to reject H would 


o(ii) 
imply that the In A(x,y).plane is tilted with respect to 
the x-y plane, but no "warping" is evidenced. 

The above development relies heavily on asymptotic 
assumptions. Small sample problems will be much more 
difficult to analyze. Any point in the above procedure w! .ch 
lead to rejection of any hypothesis would require the ana. “sis 


to proceed with the estimation of the non-zero parameters 


itis stem SUbJECT of the next chaprver- 


Di 


A "GIA 





V. ESTIMATION OF PARAMETERS 


It is desired to formulate a method for estimating the 
parameters oa, B, y and 6 of the non-homogeneous planar 
Poisson model given in IV-A where it has been established 
that a non-homogeneous process is evidenced by the data. 

Taking the logarithm of the conditional likelihood 


Bımecion (37) results in 
In L(n) = Inn! + Bix, + yly, + $Ex,y, +n1nA, (38) 


where A = A(X*,Y#)/&. Point estimation of (a, B, y, 6) 

by the method of maximum likelihood uses the conditional 
likelihood function (38) to develop the estimates. See 
Section II-D for comments regarding use of maximum likelihood 


ii Units application. The solution to the set ol simultaneous 


equations 
y* yk 
LX, -= à y J wu expíBu + yv + Suv) dudv = 0 
0 0 
y* xs 
DES —;I Í v expí8u + yv + Suv} dudv = 0 (39) 
a A 0 0 
y* x* 
LX.y. - T f uv exp{Bu + yv + Suv} dudv = 0, 
P A 0 0 ) 


A 


ratas ble, provides the polnt estimetes B, y and ô. 


Note that this approach neglects the homogeneous term 
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Meine the estimation of the parameters giving rise to 
non-homogeneity. The neglected parameter may be estimated 
last. 

In order for the solution (8, v. 5) to equations (39) 
to describe a relative maximum to 1n L|n, it is necessary 
ema sufificlent that the matrix of second partial derivatives 
(2) be negative definite, see Frisch [1966, p. 120]. In 
examining this matrix in the case of (38), it is helpful 


to define S(u,v) = exp {ßu + yv + Suv). Then the function 


s(u,v) = S(u,v) 


A 


has the properties: 
nv) > 0 
Y* Xx* 
(b) J Su) dudy e] 
0 0 


(CS O ES CONTANDO <A]. 


Hence s(u,v) is a probability density function [Gnedenko, 
"062, pre 1/1]. 
Hence the matrix 2 can be shown to have diagonal 
elements such as 
Y* x* Y* x* 


= - nl S Ji u*s(u,v) dudv - ( f jo us y) ELM 
0 0 0 0 


Q 
| 


IL 
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Continuing, the result is (where W is defined to be the 


W= UV) 
-n Var U -n Cov (U,V) -n Cov (U,W) 
r= ne (U V) -n Var V -n Cov (V,W) 
-n Cov (U,W) -n Cov (V,W) -n Var W 


Mo aS revealed to be a covariance matri Note that. the 
CORON Tor a relative maximum, Imern Dmegstivegerimite, 
is independent of the realizations. 

Now E = -n?r where E is the usual variance-covariance 
matrix Lor a tri-varíate distribution. But E SOS valve 
semi-definite [Gnedenko, 1962, p. 212], hence T is negative 
semi-definite. That each of the principal minors has 
mn zerordeterminants remains tobDeRSNOwWN. 

By the expressions given in Gnedenko [1966, p. 212], 


ag 


mise covariance matrix 2 Cam be Seeneroe Dema Hankel matrix 
[Gantmacher, 1, 1959, p. 338]. Hence if the rows of I ar 
ime ar ly independent, then the determinant of I PEE BUT 
also Var U > 0 since U is a random variable and Var U Var V - 
Cov^(U,V) » 0 since the case of line discontinuities has 

been excluded (i.e., U cannot be a linear function of V). 

By the same reasoning, W is linearly independent of U and 
presence all principal minoOmsware jgreaver than zero, hence 

z is Positive definite henc: meS mmee awe definite. 


MSE y, à) provides at leastre relative maximum to 


in Lin. 
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If it were possible to determine that (8, Y, 5) provides 
a global maximum to ln L|n in the region of interest, then 
@enclusions as to uniqueness of the estimator could be 
ieew.  Uniortunately, global “extrema are difficult to 
establish. Since the method of estimation used was maximum 
tikelihood, the estimates are consistent. Questions of 
biasedness are unresolved. 

In order to solve the system of equations (39), it is 
meeessary to determine initial values Tor che parameters as 
a Starting point for an iterative procedure. The partial 
differentiation of 1nL (35) with respect to the parameters 
and setting these partials equal to zero results, after 


some algebraic manipulation, in 


H Np = 0 
Qt yY, (Bt6Y)X ax 
Y om Cmte -1) _ ef” - 1,5 _ 
Ex; te - al B+8Y a 
BX, (ytdX)Y oa 
8 e” -e”“(e A . 
oe ee ee o o 


(BARS e? By .BX+yY+6XY _yY 8X. 
buy, * (772A * SOY) 83). 859 a de 


~ ale u - 1] - exte Y* - 1] - ES - 1]) 20 


If it is assumed that the sum BX * yY + 6XY is small 
(near zero) as well as the individual terms in the summation 


Being small, then tcthezexzprenen esse make zapproximated by 
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exp{x} 3 1 + x, x near zero. Using the first equation 
in system (40) to give the value for A(X*,Y*), i.e. 
EXC Y*) - n, and the linear approximation in the remaining 


terms gives the abbreviated system: 





e Y 
Lx. T sn = 0 
Ly, + i n = 0 (41) 
RES 

BY+S (1+0) [28 By 2 

LX Y + = n E [ot B56X* Y TROY 
E > 
-6 XY * egyxY] = 0 


Imersotucton to (41) provides the initial estimates for the 

parameters. These estimates can then Ce used in (39) or 

(40) to search for sequentially closer and closer 

approximations in a mathematical programming approach. 
Fowlowving the deternanetiemnsorzsche este mates 8, Y 

and 5, E can be determined Tromzche solution to the first 

equation in the set (40). 


The determination of confidence intervals and levels 


Sere sienaiieance is not considered. 
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VI. CONCLUSIONS 


The procedures in Chapters III-B, Iy and V are dependent 
on the particular choice of parameter form; however, with 
different forms the concept of a non-homogeneous planar 
Poisson process may be used to describe a wide variety of 
"randomly" occurring phenomena. The choice of parameters 
which may be used is limited only by assumption I, i.e. 
positivity. One advantage of the method discussed herein 
over previously proposed schemes is the fact that the ' 
specific form used admits the possibility of a ridge or 
line of maximum density to be mathematically specified 
and estimated. 

Also there is an attempt to describe the underlying 
Process Chat caused the points k oseo eadar where they did, 

Els Opposed to using, for instameesw cme are within which 
the most events were observed as the point estimate for 
the direction of maximum increase. 

Further efforts in this area include a generalization 
Mato tour dimensions (X,¥,2,0)eammenrder thal zoological 
as well as botanical densities may be studied. Of especial 
interest is the estimation of densities of aquatic life 
and how the observed density fluctuates with season and 
with changes in environment. The latter problem seems of 
prime importance in evaluating the effects of anti-pollution 
programs on the tluid systems in which plants and animals 


exist. 
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Another problem which is closely related to the above 
1S that of imperfect Sampling and how the estimates are 
biased by sampling techniques. 

Chapters III, IV and V may be redefined in terms of 
data ‘gathered within a circle about some fixed point, 
especially with consideration of the relative efficiency 


of this data form referred to by Matern [1960]. 
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APPENDIX A: THE BIVARIATEADULIBORM DISTRIBUTION 


2 of area A and the fact that the 


Given a region R in E 
probability of occurrence of an event in any sub-region R 
of area A, within R is simply ASA, a bivariate uniform 


distribution is described. For definiteness assume the 


meson R is rectangular, so A X*Y*. Now 


Prob (Xx Yk y =-= 
s x*v* 


A  EENOD (YT y) 


OS x < X5. 0 < y < Y“, in which ease icf) is apparent 
phat the coordinate axes define independently chosen 
univariate random variables. 


Also, Che density function i-mi diately 
f(x,y) » 1/X*Y* UEM oc OS 


From the density function the joint density for n indepencent 


bivariate uniform random variables is 


£( (x,y) 5--+,0x,y), 5m) = 1/ Oey 


where (x,y), denotes the i'*h pair of random variables 
selected. Now n pairs of random variables, or more simply 


n points in the plane, can only be ordered (without 
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meplacenent ) lan! ways. independemnvuol cule Ordening process 


chescne Hence, the Joint dense ion for n ordered 


varate uniform random variables ts 
£(,3) 5s Guy) (uj,n) 2 nt/QUY*)" 


where (AY) (4) is the i'P point selected in the ordering 
Neme zed, 

AS a specific example, consider the n points CO De 
Tope ARE V I C Nre Spect vosimeresspembEU cEDUOECeoTvhe y= 


component. Then 


S A and (xY) tk) = ET 


If the x-eomponents are also Ordered, then the set of 
points P = UA) (4323 isd un) defines n? points, 
of which n are known to.be "occupied," that is, to descri e 
an event. For S there exists some j such thai Y (3) 
gives the y-coordinate value on nene ment which gave rit: 
o X(1)* Sar y. for X(2) there are 

now n-1 3's remaining, one of which must correspond to the 
event giving rise to X(2)* Continuing to Xn)» phere can 
eme] left to be asseociavedmiaiunm the last x-value. 
Thus there are n! combinations of (X,Y) (4) 5 Mac] m 
each having density of 1/(X*Y*)% and so the ordered bivariate 


uniform density is established as that stated above. 
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ABSTRACT 





The occurrence of independent events at random imme plane. 
1.e. the formation of a planar Dou paces Si. S discussed. Both 
homogeneous and E processes are considerec. A 
specific functional form for the parameter in a nonhomogeneous 
planar Poisson process is used to illustrate the development of 
test and parameter estimation techniques. The problem finds 
application in the desription of biological phenomena as well as 


in search and detection problems. 
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